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ABSTRACT 



Aims. To constrain the ionization fraction in protoplanetary disks, we present new high-sensitivity interferometric observations of N2H + in 
q \ three disks surrounding DM Tau, LkCa 15, and MWC 480. 

■ Methods. We used the IRAM PdBI array to observe the N 2 H + J= 1-0 line and applied a ^-minimization technique to estimate corresponding 

' column densities. These values are compared, together with HCO + column densities, to results of a steady-state disk model with a vertical 
temperature gradient coupled to gas-grain chemistry. 
Q ■ Results. We report two N 2 H + detections for LkCa 15 and DM Tau at the 5cr level and an upper limit for MWC 480. The column density 
S—j derived from the data for LkCa 15 is much lower than previously reported. The [N 2 H + /HCO + ] ratio is on the order of 0.02-0.03. So far, HCO + 
C/3 , remains the most abundant observed molecular ion in disks. 

Conclusions. All the observed values generally agree with the modelled column densities of disks at an evolutionary stage of a few million 
years (within the uncertainty limits), but the radial distribution of the molecules is not reproduced well. The low inferred concentration of N 2 H + 
• ■ in three disks around low-mass and intermediate-mass young stars implies that this ion is not a sensitive tracer of the overall disk ionization 
' fraction. 

Key words. Stars: circumstellar matter - planetary systems: protoplanetary disks - individual: LkCa 15, MWC 480, DM Tau - Stellar 
PropertiesRadio-lines: stars 



1. Introduction 

A comprehensive understanding of the chemical composition 
and evolution of protoplanetary disks is a necessary prerequi- 
site to set up the physical and chemical conditions in which 
planet formation should occur. Apart from CO and its isotopo- 
logues, and occasionally HCO + , the molec ular content of pro- 
toplanetary disks remains poorly known. iDutrev etal] dl997l) 
first reported the discovery of a number of simple molecules 
(HCN, CN, C 2 H, H 2 CO, and CS) in the circumstellar disk 
of DM Tau and the circumbinary disk of GG Tau, using the 
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IRAM 30-m telescope, while Kastner et alj d 1997b observed 
HCN, HCO + , and CO in TW Hydra. Single-dish observa- 
tions of molecular lines have also been re ported for LkCa 15 
(Ivan Zadelhoff et alj|200lb iThi et alj|2004b and for AB Auriga 



Semenov et all (2005 



by 

More recently, millimeter arrays have started to search 
for molecular fines in the protop lanetary disks of LkCa 15 
dAikawa etal J 120031: bi et alj 12003b and TW Hydra dOi etal ' 



20041) . However, high-angular resolution observations are ex- 



tremely scarce and often noisy. Moreover, the simultaneous de- 
tection of the continuum emission from the dust disk can bias 
the determination of the line intensity if not properly taken into 
account. 

Since molecular line data are limited in sensitivity and in 
spatial resolution, little is known about the spatial distribution 
of molecules and their abundances in disks. As a consequence, 
comparison with existing chemical models are hampered and 
often based on global data (e.g. the total number of molecules 
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or, at best, integrated line profiles). With these integrated quan- 
tities, it is difficult to test the validity of sophisticated models, 
which include thermal structure and even time dependency. 

Using the IRAM PdBI array, we started in 2004 the CID 
project ("Chemistry In Disks"), a molecular survey of well- 
known disks such as those surrounding LkCa 15, MWC 480, 
and DM Tau. Our main goal was to map sufficiently strong 
molecular lines already detected in several disks with enough 
sensitivity to derive molecular abundance variations versus ra- 
dius. Among the molecules we searched for, molecular ions 
such as HCO + or N2H + are of great interest because they 
are potential tracers of the ionization structure of the disks. 
Note that the latter molecular ion is usually used as a high- 
density tracer in dense molecular cloud cores as it is supposed 
to b arely freeze out on to dust grains even at low temperatures 
(e.g. lCrapsi et al.ll2005l) . 

In this first paper of a series we report the observations 
of N 2 H + in three objects: DM Tau, LkCa 15, and MWC 480. 
Using the ^-minimization technique in the UV p lane as de 



scribed in Guilloteau & Dutrey (1998; see also IPietu et al 



20071 for details), we estimate the column density of N2FU in 
the outer parts of these disks and the corresponding uncertain- 
ties. We compare these observed values with the N2H + column 
densities computed with realistic 2-D steady-state disk mod- 
els and a gas-grain chemistry with surface reactions. We also 
briefly discuss the observed propert ies of HCO + in these disks 
(using results from Piet u et al . 2007) and deduce an estimate of 
the ionization fraction. 

In Sect. 2, we describe the observations and the sample of 
stars. The results are given in Sect. 3. Sect. 4 is dedicated to the 
corresponding chemical modelling. We discuss their implica- 
tions in Sect. 5, followed by conclusions. 

2. Sample of stars and observations 

2.1. Sample of stars 

Table[T]gives the coordinates and the physical properties of the 
observed stars. Our stars have spectral types ranging from Ml 
to A4. We chose these objects because they have large dust and 
molecular disks. Moreover, these stars are located in a hole or 
at an edge of their parent molecular cloud, in a region devoid 
of CO emission. 

2.2. Observations 

Observations of N2FU were carried out in November 2004 and 
in April and September 2005. The observing time was shared 
in each track between LkCa 15, MWC 480, and DM Tau. The 
D configuration was used with 6 antennas. The sources were 
observed in time-sharing mode during 3 transits, leading to 
a total on source integration time of about 6 hours for each 
object. At 3 mm, the tuning was single-side-band (SSB). The 
backend was a correlator with three bands of 20 MHz (spec- 
tral resolution 0.125 kms 1 ) centered on the N 2 H + lines, one 
band of 20 MHz centered on the CS J=5-4 line, and 2 bands 
of 320 MHz for the 1.3mm and 3mm continuum, respectively. 
The phase and flux calibrators were 0415+379 and 0528+134. 




MWC 480 




u 



Fig. 1. Integrated spectrum of N 2 H + J=l->0 for DM Tau 
(top), LkCa 15 (middle), and MWC 480 (bottom). The thick 
curve represents the best model deduced from the^- 2 minimiza- 
tion procedure. The location and relative intensities of the hy- 
perfine components are indicated for reference. The spectral 
resolution is 0.25 km s _I for DM Tau and LkCa 15 and is 0.5 
km s _1 for MWC 480. For DM Tau and LkCa 15, the lines are 
detected at the 5cr level, but are difficult to identify visually 
because of the seven hyperfine components. 



Observations of HCO + are reported in IPietu et al.1 42007). At 
3 mm, the primary beam of the 15-m antennas is on the order 
of 56", or ~ 8000 AU at the Taurus distance. 

For both data sets, we took great care with the relative flux 
calibration by using not only MWC 349 as a common refer- 
ence (see the IRAM flux report, number 11) but also by always 
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Source 


Right Ascension 


Declination 


Spect.Type 


Effective Temp. 


Stellar Lum. 


Stellar Mass 


Age 


UV flux 




(",',"), (J2000.0) 


(\"Y) (J2000.0) 




(K) 


(L ) 


(M Q ) 


(Myr) 


(Go) 


LkCa 15 


04:39:17.78 


22:21:03.34 


K5 


4350 


0.74 


1.01 +0.02 


3-5 


2550 


DMTau 


04:33:48.73 


18:10:09.89 


Ml 


3720 


0.25 


0.53 ± 0.03 


5 


410 


MWC 480 


04:58:46.26 


29:50:36.87 


A4 


8460 


11.5 


1.83 ±0.05 


7 


8500 



Coordinates J 2000.0 deduced from the fit of the 1.3mm continuum map of th e PdBI performed inlPietu et al. I ( 120071) . Col. 3,4, 5, and 7 are th e 
spectral type, effective temperature, the stellar luminosity a nd age as given in ISimon et al.l J2000T) . Masses are as given bv|Pietuetal.N2007r) . 
The stellar UV fluxes that are given in Col. 8 in un its of the lDrainel dl978h interstellar UV field are taken from lBergin et alJ J2004h (LkCa 15 
and DM Tau) or computed from the lKurucj dl993h ATLAS9 of stellar spectra (MWC 480). They are given for a distance of 100 AU from the 
star. 




-10 10 
R.A. Offset (") 

Fig. 2. Weighted intensity maps (see text for details) of N 2 H + J=l->0 towards DM Tau, LkCa 15, and MWC 480. The inter- 
ferometric beam is 7.0" x 4.6", with PA. 75° or 1000 AU x 650 AU at the Taurus distance. Contours are approximately 1.5 
cr. 



having our own internal reference in the observations. For this 
purpose we applied two methods: i) we used MWC 480 as a 
reference because its continuum emission is reasonably com- 
pact and bright; ii) we also checked the consistency of the flux 
calibration using the spectral index of DM Tau, as described in 
Dartois et al.l d2003l) . We cross-checked both methods and thus 
obtained a reliable relative flux calibration from one frequency 
to another. Using our method, we estimated the relative accu- 
racy of our flux density scale to be significantly better than 10% 
at 3 mm. 

3. Results 

Integrated spectra (over a 6" x6" region) are displayed in Fig[TJ 
superimposed with the best-fit disk model, as derived from the 
disk modelling described below. 

The N 2 H + J=l— >0 emission is very weak, and determina- 
tion of the integrated area is complicated by the hyperfine struc- 
ture. Furthermore, due to the Keplerian velocity gradient, some 
components overlap in velocity. Half of the signal is in the 
satellite components: computing the integrated line flux would 
require us to include these low S/N data using twice more band- 
width than for the main group. This would degrade the overall 
S/N ratio. To better illustrate the detection, we instead applied 
a filter in which each channel is weighted by the expected line 
intensity (obtained from the best fit), and the sum over channels 
was normalized by the sum of weights, a classical optimal fil- 



tering. This method confirms the detection but does not provide 
a quantitative estimate, which is given by the minimization (see 
below). The resulting weighted intensity maps are displayed in 
Fig|2] Weak but significant detections are visible at the posi- 
tion of DM Tau and LkCa 15, but no signal is detected towards 
MWC 480. 

The column density of N 2 H + is not a directly observable 
quantity, but it must be deduced from the map intensities, tak- 
ing the line formation mechanism into account. In the analysis 
of disk line data the surface density is often assumed to be a 
power law, and then the power law index is fixed during the 
analysis. In this paper we follow a different approach by ap- 
plying a minimization technique that results both in the N 2 H + 
column density and the power law index. 

Since the analysis of the N 2 H + J=l— >0 line is complicated 
by its hyperfine structure, we decided to estimate the N 2 H + col- 
umn density by directly compari ng the interferometric data t o 
molecular disk models. Following Gui lloteau & Dutreyl(ll998l) . 
the data were analyzed inside the UV plane. For each source, 
we fit a passive Keplerian disk model that assumes power- 
law radial dependencies for the molecular rotation temper- 
ature (T r {r) = 7/100 x (ioHau) 9 ) an d the surface density 
(iV N2 H+M =iViooX( 



100AU 



100AU •* 

)~ p ). We assume hydrostatic equilib- 



rium with a scale height derived from the CO measurements 

I 1 1 1 

(see Pietu et al. 2007), and a fixed local line-width Av. The 



model also takes into account the stellar mass (Table [TJ and 
disk inclination (Table |2). The disk parameters were taken as 
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o 1 

Column density exponent 

Fig. 3. Surface density of N 2 H + at 400 AU (with error bars) 
versus radial exponent for DM Tau, the error bars are for +1 
formal cr. The units are specified on the left Y axis. The light 
grey curve is the y 1 (with an offset), with units on the right Y 
axis. 
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1 
Column density exponent 

Fig. 4. Surface density of N 2 H + at 300 AU (with error bars) 
versus radial exponent p for LkCa 15, the errors bars are for +1 
formal cr. The units are specified on the left Y axis. The light 
grey curve is the x 1 (with an offset), with units on the right Y 
axis. 



constant during the^ 2 minimization process. Only the surface- 
density law parameters (A^ioo and p) of N 2 H + were varied dur- 
ing the minimization process. For thermalized lines such as CO 
J=l-0 and J=2-l or HCO + J=l-0, the rotation temperature is a 
direct estimate of the gas kinetic temperature. In our case, we 
take the rotation tem perature law from HCO + data (see TablefJ] 
and lPietu et al.ll2007l for details). Such a low temperature im- 
plies that the HCO + J=l-0 emission is coming from the mid- 
plane. Given the expected temperature range that is measured 
in these disks (see Tables |2][3]), the chosen values have little 
influence on the derived column density. 

The results of the fits are summarized in Table [3] and pre- 
sented in Fig. [3] for DM Tau and in Fig. [4] for LkCa 15, re- 
spectively. We have chosen to show the results at 400 AU for 
DM Tau (300 AU for LkCa 15) because it is the relevant scale 
for the IRAM data in terms of sensitivity and resolution. These 
figures indicate that we detected N 2 H + at the 4 - 5 <x level in 
both sources, with a best-fit exponent of 2 - 2.5. For MWC 480, 
the formal best fit leads to negative column densities, at the lcr 
level: we thus only mention a 3<x upper limit for this source. 

The 5cr detection in LkCa 15 corresponds to a column den- 
sity that is apparently about a factor ~ 100 below th e column 
density deduced from OVRO data by lOi et aL (2003), who re- 
port an observed value of 1.7 1 13 cm -2 . The difference can be 



Table 2. Disk Parameters 



traced down to several factors. lOi et alj d2003l) indicate a total 
line flux of 0.3 Jy km/s (from their Tabl e 1), whi l e our best fit 
line flux is 0.15 Jy km/s. Furthermore. lOi et alJ (120031) use a 
constant temperature of 30 K, and give the beam-averaged col- 
umn density (p = 0), which corresponds to a source radius of ~ 
280 AU. Under the same conditions, we deriv e a column den- 
sity of 4.5 + 0.9 10 12 cm" 2 . Assuming that the bi et all d2003l) 
flux corresponds to the main group of hyperfine components 
resolves the discrepancy and shows we have detected 4 times 
less N 2 H + molecules than claimed by lOi et all (12003 ) if we 
take the same conditions for the analysis. Part of the difference 
may be due to the treatment of the continuum: when comput- 



Source 


i 




^100 


Acc.Rate 


Disk Mass 




o 


(AU) 


(K) 


(M yr-') 


(M ) 


DM Tau 


-32 


800 


15, q~ 0.5 


2.10- 9 


0.05 


LkCa 15 


52 


680 


15, q~ 0.5 


io- 8 


0.03 


MWC 480 


38 


695^ 


30, q~ 0.5 


io- 8 


0.03 



This value is taken from lThi et aD J2004h . 



Col. 2 Inclinations are taken from ISimon et al.1 feOOd) . Col.3 & 4 Outer 
radius and temperature a re taken from the 13 CO interferometric anal- 
ysis of lPietu et alJ l l2007h . Note that the temperatures given here corre- 
spond to the observed (CO data) gas kinetic temperature from the cold 
mid-plane region and is higher in the disk intermediate layer. They are 
in good agreement with predictions from the model by D' Alessio and 
collaborators. Col.5 A ccretion rates are taken as a generic value from 
ID' Alessio etafl dl999h for MWC 480 and LkC a 15. The accretion rate 
for D M Tau is taken fromlCalvet et a l. (2005). C ol.6 Di sk masses are 
fro mlDutrev et alj J 19971) for DM Tau. iThi .et ail J2004I) for LkCa 15 
and lMannings et al.Ul997l) for MWC 480. 



ing the integrat ed emission, we subtracted the continuum while 



Oi et al.l d2003l) included it with the justification that the lines 
are optically thick. The assumption that the lines are optically 
thick clearly does not a gree with intensitie s detected so far i n 



disks, as HCO + (1-0) dSemenov et al.ll2 005; Pie tu et al.ll2007l) . 
C 18 O(l-0) and 13 CO(T-0) dDartois et al.ll2003l) . and NoH+ri-OI 



are optically thin. Furthermore, even for an optically thick line 
the continuum should be subtracted because at a given veloc- 
ity in the disk, the line emission covers only a tiny fraction of 
the entire continuum, which is at most t/V/VKepier where dV is 
the local line-widt h and Vxepier the projected rotation velocity 
at the disk edge ( Beckwith & Sargentf 19931: iGuilloteau et al 



2006). Thus the integrated area of the line will be significantly 



overestimated if the continuum is not subtracted. 
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Table 3. Column densities of HCO + and N 2 H + and abundance 4.1. Model description 

ratios 



Source 


MWC 480 


LkCa 15 


DMTau 


T 


20 


19 


13 


q 


0.4 


0.38 


0.15 


N(HCO + ) 


5.0±0.610 12 


8.0 ±0.5 10 12 


6.5 ± 1.4 10 12 


P 


2.1 ±0.5 


2.5 ±0.3 


2.7 ± 1.0 


R 


250 AU 


300 AU 


400 AU 


JV(N 2 H + ) 


<310" 


2.3 ±0.5 10" 


1.1 ±0.25 10" 


P 


2 


2.5 


2.5 


R 


250 AU 


300 AU 


400 AU 


[N 2 H + ]/[HCO + ] 


<0.06 


-0.03 


-0.02 



T indicates the temperature (in K) at 100 AU and q the t emperature 
expon ent used to derive the column densities, taken from Pie tu et al.l 
(2007). For DM Tau and LkCa 15, in agreement with the figures, we 
present the results for p = 2.5. Columns de nsities are in cm" 2 . Values 
for HCO + are taken from lPietu etafl d2007l) . 



The l arge discrepancy between our analysis and that of 



Oi et al.l (120031) clearly demonstrates the necessity of taking 
proper care of the column density distribution. When compar- 
ing molecules with different apparent extents or observed with 
different interferometric beams, this is fundamental in order to 
avoid biases due to sensitivity limitations. In our case, we found 
p = 2.2 + 0.8 for DM Tau and p = 1.6 ± 0.6 for LkCa 15. These 
values are deduced from the best N2H + models and are in com- 
plete agreement with the analysis performed for the most sensi- 
tive tracer HCO + . Most importantly, we derive [N 2 H + ]/[HCO + ] 
under the assumption that this ratio is constant with radius. This 
is justified both by the above agreement between the p values 
and by the chemical models presented in Sect.4. 

This result should serve as a warning that "observed" col- 
umn densities either (i) come with assumptions for the surface 
density profile or (ii) are based on a fitting procedure within 
physical models. In our work, the terms "observed column den- 
sities" refer to values derived using the second method. The 
expression "modelled column densities" refers to column den- 
sities predicted by a chemical model, such as the one described 
in the next section. 



4. Chemical modelling 

Since our interferometric millimeter data are sensitivity lim- 
ited, we have chosen to focus on the more prominent chemical 
effects. In this paper, our main goal is to reproduce the over- 
all distribution of the observed molecular ions. Accordingly, 
the chemical model we use does not t ake into account grain 
growth and sedimentation. Similarly, ISemenov et al ] d2006l) 
have shown that turbulent transport can modify the column 
densities and abundance distributions of some molecules with 
respect to the usual static chemical models in a number of 
cases. However, we do not include vertical/radial mixing here. 



In order to confront the observed values with theoreti- 
cal predictions, we used an updat ed gas-grain chemical 



model with surface reactions from Se menov et al.l (12005). 



Brie fly, it is base d on the UMIST'95 gas-phase network 
( Mil lar et al. 1997 ) supplemented with surface chemistr y on 
uniform O.lyum silicate grai ns from lHasegawa et alj (11992) and 
Hasegawa & Herbstl (119931) and a small set of deuterium chem- 
istry reactions (E.A. Bergin, private communication). Using 
measurements published during the last decade, we updated a 
few tens of important reaction rates in the UMIST'95 database 
(IVasvunin et a 1. 2006). Overall, our chemical network consists 
of 560 species made of 14 elements and about 5400 reactions. 

The primary ionization sources in this model include stel- 
lar X-ray and UV radiation, interstellar UV radiation, and cos- 
mic ray particles, and the decay of short-lived radionuclides. 
The adopted intensity of the un-attenuated stellar UV flux 
at 100 AU from the central stars are given in the last col- 
umn jjfTableffl We used the standard interstellar UV field 
of iDraind d 19781) and re-scaled it to match the correspond- 
ing stellar luminosities of TTauri stars according to the current 
limited knowledge of their UV spectra. This approach should 
be improved as so on as better UV data become available (see 
Bergin et al.l l2006). Note that the UV penetration is computed 
using a 1-D plane-parallel approach and thus the UV inten- 
sity in a given disk location is very likely underestimated in 
com parison with the full 2-D radiative transfer simulations; 



van Z adelhoff et al. (2003). The intensity and penetration 



of the stellar X-ray radiation are m odelled using observa- 
tional resul ts of Gl assgold et al. (2005) and Monte-Carlo sim- 
ulations of iGlassgold et al. ( 1997a b). The corresponding un- 
attenuated ionization rates at 100 AU are equal to 4.5 10~ 14 s _1 



for DM Tau, 
MWC 480. 



2 10-" s" J for LkCa 15, and 2 lfT 



for 



Using the 



resul ts of laboratory measurements by 
Bisschop et al.1 (120061 their "mixed-ice" case), we assumed 
that desorption energies of CO and N 2 are equal to 930 K. For 
all sources, we used a flared 1 + 1D di sk model with a vertica l 
temperature gradient similar to that of ID' Alessio et al. I dl999h . 
This disk model has the - 30 K temperature surface at - 3 
pressure scale heights (the CO surface where r - 1 for the CO 
J=2-l transition) at a di stance of 100 AU, in agre ement with 



the re sults obtained by iDartois eta l. (2003) and IPietu et al 



(120071) . Other disk parameters are summarized in Table[2] 



With this physical and chemical model in hand and TMC1- 
like initial abundances, the chemical evolution is simulated 
over the 5 Myr (DM Tau and LkCa 15) and 7 Myr (MWC 480) 
evolutionary time. We repeated the modelling of the disk chem- 
ical evolution 20 times by using the same chemical network and 
varyin g its rate coeffic i ents w ithin their uncertaint y limits, sim- 
ilar to lvasyunin etalJ d2004l) and lWakelam et al.l d2005l) . This 
allows us to get first-order estimates of the intrinsic uncertain- 
ties of the theoretical abundances. 
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DMTau LkCal5 MWC480 




R, AU 



Fig. 5. From top to the bottom: model-averaged relative abundance distributions (in respect to the total number of hydrogen 
nuclei) of important ions: H + , C + , H+, HCO\ N 2 H + , and CO in the DM Tau (left), LkCa 15 (middle), and MWC 480 (right 
panel) disks after 5 and 7 Myr of the evolution. The Y-axis is the disk vertical extent normalized by one pressure scale height 
(different for all three objects). 



4.2. Ionization in disks and global evolution of N 2 H + 
and HCO + 



Y-axis is expressed in units of the pressure scale height, H r , as 
defined in lDartois et al.l ([2003). 



More generally, in the upper surface layer (k, 2 scale 
In Fig. [5] we present the model-averaged disk abundances of heights) the dominant ion is H + , with an abundance up to 
CO, several dominant ions, and the ionization fraction after a ~ 10~ 3 (with respect to the total amount of hydrogen nuclei), 
few Myr of evolution for all three objects without mixing. The Somewhat lower, at ^ 1.5 scale heights, the disk ionization de- 
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gree is determined by ionized carbon, which is also the most 
abundant ion in all three disks ([C + ] ^ 1(T 4 , N(C + ) ~ 2V(e _ ) ~ 
3 10 16 cm 2 around ~ 200-500 AU). The evolution of both 
ions is governed entirely by fast photo-reactions due to ener- 
getic stellar UV field and X-rays and interstellar UV radia- 
tion. In contrast, the less abundant H^ ion mostly traces the 
fractional ionization closer to the disk midplane, [HJ] ~ [e~] 
$ 10~ 9 , where the cosmic ray particles remain the only ioniza- 
tion source. Unfortunately, all these important ions are not eas- 
ily detectable in protoplanetary di sks (however, see lKamp et al 
120031 ICeccarelli & Dominikll2005l) . 

The observationally most important ion in the outer disks is 
predicted to be HCO + (typical abundance value is % 10~ 9 ). It is 
located in the warm molecular layer that is partly shielded from 
direct ionizing radiation from the central star and the ISM, at 
~ 1 pressure scale height. Note that the corresponding HCO + 
layer is thinner and located deeper towards the midplane in 
the MWC 480 disk in comparison with DM Tau and LkCa 15. 
Such a chemical effect is caused primarily by the strong radia- 
tion field of MWC 480, which is capable of dissociating CO - 
the key molecule for the HCO + evolution - more deeply inside 
the disk, as seen in the figure. The more intense stellar heating 
also leads to a warmer interior of the MWC 480 disk compared 
to the cases of DM Tau and LkCa 15, which in turn results 
in efficient CO desorption in the MWC 480 disk midplane at 
R £ 200 AU and thus in a reservoir of abundant HCO + there. 
Note that, in all three cases, the upper boundary of the abun- 
dant HCO + layer coincides with the lower boundary of the C + 
layer that marks the disk height where penetration of the ioniz- 
ing radiation essentially stops. This happens at ~ 1 — 1.5 scale 
heights. 

In contrast to HCO + , the next observationally important ion 
N2H + is concentrated in a very thin layer at the lower bot- 
tom of the intermediate layer (~ 0.5 - 1 scale heights) and in 
the disk midplane. Unlike for other considered species, N2H + 
abundances are more spatially peaked toward outer disk re- 
gions. It is also interesting that the N2H + vertical layer is in 
general thinner in the MWC 480 disk compared to the disks 
around LkCa 15 and DM Tau, similar to the HCO + distribu- 
tion. The main reason for the chemical differentiation between 
these two ions is the fact that N2FU gets destroyed easily in 
proton transfer reactions wi th CO, e.g. N 2 H + + CO - » HCO + 
+ N 2 (k = 
Thus, 



NH3 by a two-step process of surface hydrogenation: 



- 10 cm 3 s- 1 



Anicich & Huntressfl 986) . 



8.810" 

N2H + is less abundant in the disk region where the 



CO concentration is high, that is, in the CO molecular layer 
above about one pressure scale height. On the other hand, it 
also cannot be abundant in the disk surface layer exposed to the 
energetic radiation from the star, which destroys all polyatomic 
species within ~ 10 4 years. Finally, it cannot reach high levels 
of concentration in the outer disk midplane (R > 400 AU in 
DM Tau and LkCa 15, and R > 200 AU in MWC 480) because 
the evolution of N2H + in the entire midplane region is governed 
by a cycle between ion-molecule reactions of N2 with the most 
abundant ion, H^ (5= 10~ 9 ), and dissociative recombination of 
N2H + with electrons and charged grains. Consequently, at very 
low temperatures in the outer midplane a volatile product of the 
N2H + dissociation, namely NH (desorption energy is 600 K), is 
frozen out on grain surfaces and eventually converted to solid 



NH + H grains NH 2 + H grains NH 3 . 

Thus, after about 1 Myr of evolution all nitrogen is locked in 
solid ammonia that has an abundance of about 10~ 5 . In con- 
trast, in the warmer disk midplane region that is located closer 
to the central star, NH is not totally frozen out and thus the 
surface formation of ammonia is not that efficient. Therefore, 
about 10% of all nitrogen is locked in solid N2 (~ 2 10~ 6 ), but a 
minor fraction of molecular nitrogen remains in the gas phase 
(~ 10~ 8 ) and reacts with H+ , forming N 2 H + . 

Note that there is a thin layer devoid of N2H + at about 0.7 
scale height in the disk of DM Tau at 100 AU and, to a lesser 
extent, LkCa 15. In contrast to the Herbig Ae star MWC 480, 
these two Sun -like T Tau stars are s trong emitters of hard X- 
ray radiation dGlassgold et al. I l2005l) . The stellar X-ray radia- 
tion cannot propagate directly through the disk midplane, but 
in general has larg er penetration depth toward the disk interior 



than UV photons. Sem enov et al.l (120041) have shown that the 



attenuated X-rays can drive a peculiar chemistry leading to the 
production of complex (organic) species on dust grains. 

5. Discussion 

5.1. Observed abundances of molecular ions in disk 

HCO + is observed in many astrophysical media with a wide va- 
riety of physical condition s, from diffuse clouds (ILiszt & Lucas 



2001) up to dense cores (Casell i et al.ll2002al) . N 2 H + , which 



can be signifi cantly less abund ant and even depleted in very 
dense cores (Pagani et al.l l2005l) . appears more difficult to ob- 
serve dTafalla et al.l 2004). Moreover, as described in Sect. 4.2, 
N2H + c an be destroyed by re actions with CO to give HCO + 
(see also. lAnicich & Huntressl[l986]) . Hence, HCO + and N 2 H + 
usually do not coexist. Having abundance estimates coming 
from the same area is not so easy. This is, for example, illus- 
trated by J0rg ensenl ( 2004 ) in L483-mm where he finds that the 
HCO + traces the large scale outflow, while N2H + traces the 
quiescent material (in the core) being accelerated by the out- 
flow (see also his Fig. 13). Observations of abs orption lines in 
diffuse clouds in front of bright quasars allow Liszt & Lucas 
(2001) to estimate a lower limi t on the abundance ratio of 
[N 2 H + ]/[HCO + ]< 0.002 (at 2cr). lAikawa eTaD (120051) predict 
abundance ratios of [N 2 H + ]/[HCO + ] within the range 0.02-0.1 
for physical conditions relevant for pre-stellar collapsing dense 
cores. These chemical models correspond to abundances of 
N2H + and HCO + , in agreement with observations of starles s 
cores such as L1544 dTafalla et al.ll2004t ICaselli et al.ll2002bl) . 
To first order, we get an averaged ratio [N2H + ]/[HCO + ] on the 
order of ~ 0.02 - 0.03 (see Table |3). The abundance ratios we 
observe in the TTauri disks are therefore similar to those ob- 
served and modelled in dense cores of cold molecular clouds. 



5.2. Reliability of chemical model 

Our non-mixing chemical model (described in Sect. 4.2) cou- 
pled to a 1 + 1D steady-state disk model with small sub-micron 
size dust grains cannot reproduce the complexity of the disk 
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Fig. 6. Top: Expected column density of N2H + versus radius (black dots), compared to the observed column densities (curve 
and en-orbars) in MWC 480 (left), DM Tau (middle), and LkCa 15 (right). Bottom: as above but for HCO + . The disk age is 5 
Myr for DM Tau and LkCa 15 and 7 Myr for MWC 480. The abrupt changes in predicted column densities of N2H + at some 
radii are partly due to the limited resolution of the chemical model (see Fig . O ■ For HCO + , the error bars include the (measured) 
uncertainty on the excitation temperature for DM Tau, while this excitation temperature was imposed for LkCa 15 and MWC 480 



structure. Surprisingly, we are able to obtain reasonable global 
agreements with the observed column densities of N2H 4 " and 
HCO + around 400 AU in all three disks (see Fig. [6). But we 
fail to reproduce in detail the observed slopes p. The modelled 
N2H + column density is in general too high in the outer disk 
(R > 400 AU), while the HCO + column density appears too 
low in the inner part of the disks. However, in the model as 
in the observations, the [HCO + ]/[N2H + ] ratio is almost con- 
stant with radius. More precisely, the value of ~ 3 10 11 crrT 2 
for the column density of N2FL in the DM Tau disk is about 
50% higher than the observed value, which is well within the 
error bars of a factor of 3 in the abundance due to the in- 
accuracy of the reaction rates in the UM IST'95 ratefile (see 
Vasvunin et al ] |2004tlWakelam et alj|2 005). The overall agree- 
ment between the modelled and observed column densities of 
HCO + in DM Tau, LkCa 15, and MWC 480 is not as good as 
for N2FL (see Fig.|6]l. The observed value of about 5 10 12 cm" 2 
that is nearly equal for the three disks is 3 times higher than 
predicted by the MWC 480 model and 3 times lower than ob- 
tained with the DM Tau and LkCa 15 disk models. Note that 
the computed HCO + column densities also suffer from intrin- 
sic uncertainties caused by the reaction rate uncertainties, about 
a factor of 3-4. 

Since the chemistry of HCO + is sensitive to the grain 
growth, the presence of large grains, as usually observed in 
disks with dust spectral index /3 - 1 - 0.5, may also account 
for the discrepancy between the results of the applied chemical 
models and observational data, especially for HCO + , but due to 



the limitations of the current chemical prescription this cannot 
be easily checked. 

5.3. Determining the disk ionization fraction 

As shown in the previous section, the ionization fraction in dif- 
ferent disk regions is dominated by distinct ions, with C + the 
most abundant one in the upper layer. Is there any possibility 
of estimating the degree of ionization in some disk regions or 
even in the entire disk? 

The HCO + abundance gives a lower limit to the electron 
fraction, [e~] > [HCO + ]. The H2 column density is uncer- 
tain: taking [ 13 CO ] = 10~ 7 and column density val ues from 
IPietu et all d2007l) leads to [e ] > 2 10~ 10 . Following loi et al 



(2003, their Eq.6), a lower limit to [e ] can also be obtained 
from N 2 H + : 



[e~] > [CO] 



[N2H+] =4.3 10-^.[CO] 



k rec (UCO + ) [HCO + ] 



[Hccr 



(1) 

where £2 = 8.8 10 -10 cm 3 .s _1 is the reaction rate coefficient of 
CO with N 2 H + (CO + N 2 H + -> HCO + + N 2 ) and fc rec (HCO + ) = 
3 10~ 7 (77300)~ a64 cm 3 s _1 is the dissociative recombination 
rate coefficient of HCO + . The numerical value is for T — 15 K. 
Thus N2FL can be a more sensitive tracer of the ionization frac- 
tion in disks than HCO + only if 



[N 2 H + ] 



1 [HCO + ] 



[HCO + ] 4.3 10- 4 [CO] 



(2) 
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Table 4. Estimated ionization degrees 



Source 


T 


n-a 


(x 


[e-] 


[e-] 




(K) 


cm" 3 


s- 1 


(HCO + ) 


(N 2 H + ) 


DMTau 


15 


5 10" 


410-' 7 


2 IO" 9 


4 10-" 


LkCa 15 


15 


5 10" 


io- 16 


3 10-' 


7 10- 9 


MWC 480 


25 


5 10 6 


1.3 10~ 17 


2 IO" 9 


310-' 



Estimates of the disk ionization degree in the intermediate layer de- 
rived from chemical equilibrium approach (EqfJJ, where the two val- 
ues (col 5 and 6) depend on whether HCO + or N 2 H + dominates the 
recombination process. 



From IPietu etail d2007t) . the [ 13 CO]/[HCO + ] ratio is 600 
After correction of the isotopologue column de nsity ratio 
around 20 when fractionation is accounted for dPietu et al 



2007), we find that, if observed with a similar S/N, N2H + can 
be a more sensitive tracer than HCO + only if [N 2 H + ]/[HCO + ] 
> 0.2. With the measured abundance ratios, ~ 0.02, HCO + thus 
remains a better tracer of the electron abundance than N 2 H + . 
Furthermore, these low abundance ratios no longer support a 
high value for the gas phase [N 2 ]/[CO] ratio. 

Another simple approach, but independent of the observa- 
tions, to obtain an estimate of the disk fraction al ionization is t o 
use the equation of chemical equilibrium (e.g.. iGammiel 1 9961) : 



[e-] = 



£rec«H 



(3) 



where f is the ionization rate (in s ), k lec is a typical re- 
combinatio n rate (cm 3 s" 1 ) , and «h is the hydrogen num- 



ber density. ISemenov et al.l (2004) have demonstrated that the 
equilibrium approach is appropriate in most disk regions. 
For the HCO + and N 2 H + ions, the recombination rates in 
our chemical network are k rec = 3 lQ- 7 (?7300) - a64 cm 3 s" 1 
and fc rec = 10- 7 (77300)~°; 5 cm 3 s -1 , respectively dMillar et al 



1997tlGeppert et al.l2004ll2005b . The steady-state values of the 
ionization degree in the warm molecular layers of the outer 
disk range are presented in Table |4] Tempe ratures were taken 
from Table|2] the density fromlDartois et alj d2003l) . and the X- 
ray ionization rates from iGlassgold et al.l (12005b . The cosmic 
ray rate is 1.3 IO -17 s _1 . The values in Table [4] are consistent 
with the results of the detailed chemical calculations shown in 
Fig. |5] (first row, one scale height and radius ~ 200-500 AU). 
These values are high enough to allow magneto-rotational in- 
stability to operate in this region of the disks (see discussion in 
Semenov et alj]2004 . 



6. Summary 

Using the IRAM array, we report new observations of N 2 H + 
J=l-0 in three protoplanetary disks surrounding two T Tauri 
stars: LkCa 15 and DM Tau and the Herbig Ae/Be star 
MWC 480. We analyze these data by also taking the HCO + 
J=l-0 data from lPietu et al.l d2007t) into account. For the three 
objects, we derive from the observations the column densities 
and abundances and compare them to model predictions ob- 
tained using a 1 + 1D chemical code. Our main conclusions are 
the following: 



- We detected N 2 H + (at 5<x level) in the disks surrounding 
LkCa 15 and DM Tau. In the disk of LkCa 15, the col- 
umn density we observed is a factor ~ 100 below previous 
claims. 

- The N 2 H + column densities derived from the observations 
imply an [N 2 H + ]/[HCO + ] ratio which is on the order of the 
ratio found in cold dense cores ~ 0.02 - 0.03. 

- HCO + remains the more abundant molecular ion in disks. 
A lower limit of 2 IO -10 for the ionization fraction can be 
obtained from the observations of HCO + and 13 CO. 

- The overall column densities of N 2 H + and HCO + and the 
[N 2 H + ]/[HCO + ] ratio are approximately reproduced by the 
chemical models, but their radial dependencies (slope p) 
are not. The latter disagreement may be linked to the lim- 
its of the chemical modelling, which does not include pro- 
cesses like grain growth or turbulent mixing. It could also 



be due to an over-simplistic disk structure, since Piet u et al 



(2006) tave shown that the dust distribution in LkCa 15 and 
MWC 480 are quite different from the simple density dis- 
tribution adopted here. 
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ABSTRACT 



v f*. ' Aims. To constrain the ionization fraction in protoplanetary disks, we present new high-sensitivity interferometric observations of N 2 H + in 
q I three disks surrounding DM Tau, LkCa 15, and MWC 480. 

■ Methods. We used the IRAM PdBI array to observe the N 2 H + J= 1-0 line and applied a ^-minimization technique to estimate corresponding 

' column densities. These values are compared, together with HCO + column densities, to results of a steady-state disk model with a vertical 
temperature gradient coupled to gas-grain chemistry. 
Q ■ Results. We report two N 2 H + detections for LkCa 15 and DM Tau at the 5cr level and an upper limit for MWC 480. The column density 
S—j derived from the data for LkCa 15 is much lower than previously reported. The [N 2 H + /HCO + ] ratio is on the order of 0.02-0.03. So far, HCO + 
C/3 , remains the most abundant observed molecular ion in disks. 

Conclusions. All the observed values generally agree with the modelled column densities of disks at an evolutionary stage of a few million 
years (within the uncertainty limits), but the radial distribution of the molecules is not reproduced well. The low inferred concentration of N 2 H + 
• ■ in three disks around low-mass and intermediate-mass young stars implies that this ion is not a sensitive tracer of the overall disk ionization 
' fraction. 

Key words. Stars: circumstellar matter - planetary systems: protoplanetary disks - individual: LkCa 15, MWC 480, DM Tau - Stellar 
PropertiesRadio-lines: stars 



1. Introduction 

A comprehensive understanding of the chemical composition 
and evolution of protoplanetary disks is a necessary prerequi- 
site to set up the physical and chemical conditions in which 
planet formation should occur. Apart from CO and its isotopo- 
logues, and occasionally HCO + , the molecular content of pro- 
toplanetary disks remains poorly known. ? first reported the 
discovery of a number of simple molecules (HCN, CN, C2H, 
H2CO, and CS) in the circumstellar disk of DM Tau and the 
circumbinary disk of GG Tau, using the IRAM 30-m telescope, 
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while ? observed HCN, HCO + , and CO in TW Hydra. Single- 
dish observations of molecular lines have also been reported 
for LkCa 15 (??) and for AB Auriga by ?. 

More recently, millimeter arrays have started to search for 
molecular lines in the protoplanetary disks of LkCa 15 (??) and 
TW Hydra (?). However, high-angular resolution observations 
are extremely scarce and often noisy. Moreover, the simultane- 
ous detection of the continuum emission from the dust disk can 
bias the determination of the line intensity if not properly taken 
into account. 

Since molecular line data are limited in sensitivity and in 
spatial resolution, little is known about the spatial distribution 
of molecules and their abundances in disks. As a consequence, 
comparison with existing chemical models are hampered and 
often based on global data (e.g. the total number of molecules 
or, at best, integrated line profiles). With these integrated quan- 
tities, it is difficult to test the validity of sophisticated models, 
which include thermal structure and even time dependency. 



